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The  results  of  temperature  measuring  experiments  conducted  for 
Rainier  shot.  Operation  Plumbbob  are  described.  The  temperature  distri- 
bution in  the  surrounding  tuff  resulting  from  the  detonation  of  an  under- 
ground nuclear  device  yielding  1.  7 kilotons  of  energy  has  been  measured. 

Data  indicate  peak  temperatures  in  the  vicinity  of  90*C  existing  in  the  cen- 
tral regions  which  drop  rapidly  to  approach  ambient  temperatures  at  dis- 
tances of  about  180  feet  from  the  original  ground  aero.  Based  on  meas- 
urements taken  in  three  holes  drilled  into  the  central  regions,  contours  of 
constant  temperature  are  constructed.  The  temperature  picture  is  consist- 
ent with  the  consideration  that  most  of  the  heat  entered  the  water  contained 
originally  in  the  tuff  and  the  unconsolidated  material  which  filled  the  cavity 
after  its  initial  collapse.  The  amount  of  energy  contained  in  the  central 
regions  in  the  form  of  tuff  heated  to  below  90*C  is  roughly  estimated  to  be 
greater  than  one-half  of  the  total  energy  release.  An  unsuccessful  attempt 
to  measure  the  temperature  rise  across  the  shock  front  of  the  pressure  wave 
produced  by  the  detonation  is  described,  and  reasons  for  failure  are  discussed. 
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CHAPTER  I 
INTRODUCTION 


1. 1 OBJECTIVES 

The  objective!  envisaged  for  this  project  were  threefold  in  nature. 

a.  It  was  hoped  that  the  temperature-time  history  of  the  tuff  surround- 
ing the  shot  cavity  would  be  helpful  in  establishing  a yield  figure  for  the  nu- 
clear device  tested. 

b.  Temperature  measurements  would  be  useful  in  documenting  the 
gross  effects  of  the  nuclear  explosion  underground. 

c.  Suggestions  had  been  made  that  the  heat  trapped  in  the  ground 
might  be  useful  for  power  production,  mining  or  petroleum  production, 
and  so  temperature  distributions  would  be  required  for  these  feasibility 
studies. 

Several  additional  developmental  programs  were  undertaken.  The 
first  involved  instrumenting  the  blast  doors  of  the  underground  access 
tunnel  with  pressure  indicators  to  determine  post-shot  pressure  condi- 
tions in  the  sealed  tunnel.  The  second  involved  the  development  and  in- 
stallation of  rapidly  responding  temperature  sensors  to  read  the  tempera- 
ture rise  across  the  shock  front  of  the  pressure  wave  generated  by  the 
detonation.  These  measurements  would  be  useful  in  establishing  the  equation 
of  state  for  the  tuff  and  in  evaluating  the  rate  of  energy  loss  in  the  medium 
as  a function  of  distance  from  the  blast. 

The  nature  of  the  program  changed  several  times  during  the  course 
of  its  execution,  but  in  general  it  was  concerned  with  pressure  and  tem- 
perature measurement.  It  is  therefore  of  interest  to  survey  briefly  the 
anticipated  pressure  and  temperature  effects. 


1.2  BACKGROUND 

The  program  was  undertaken  as  Project  26.  3 of  shot  Rainier  of 
Operation  Plumbbob  in  1957.  The  Rainier  shot  (Sept.  19*  1957)  was  an 
underground  detonation  yielding  approximately  1.  7 kilotons,  set  off  in  a 
chamber  approximately  900  ft  below  the  surface  of  a meaa  in  the  Nevada 
Proving  Grounds.  Access  to  the  chamber  was  through  a tunnel  starting  at 
the  side  of  the  mesa  about  1700  ft  from  the  chamber,  containing  several 
turns  in  the  vicinity  of  the  chamber. 

As  it  appeared  likely  that  this  would  be  the  first  of  a series  of  under- 
ground shots,  it  was  desired  that  the  instrumentation  be  as  complete  as 
possible.  At  the  same  time,  it  was  foreseen  that  difficulties  might  develop, 
and  it  was  hoped  that  these  would  be  overcome  and  not  carry  over  to  future 
underground  shots. 
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At  the  instant  of  firing  of  the  shot,  a large  quantity  of  energy  is  gen- 
erated. Since  the  radiation  is  contained  within  the  mesa,  essentially  the 
entire  amount  of  energy  is  rapidly  degraded  to  heat  and  mechanical  energy. 

It  was  believed  that  a cavity,  roughly  spherical  in  shape,  would  be  formed 
around  the  point  of  origin  by  vaporisation,  melting  and  compaction  of  the  rock. 
Inside  the  cavity  it  was  expected  that  a high  temperature,  in  the  neighbor- 
hood of  the  melting  point  of  the  rock,  would  initially  prevail.  A high  pres- 
sure was  also  expected  in  the  cavity  due  to  presence  of  gases  evolved  from 
the  rock  at  high  temperatures.  These  might  include  occluded  gases,  vo- 
latiles such  as  water,  and  other  volatile  materials  produced  in  high-tem- 
perature chemical  reactions.  A shock  wave  was  expected  to  propagate 
into  the  mesa,  with  consequent  crushing  and  deformation  of  the  rock  struc- 
ture. With  increased  distance  from  the  originating  point,  this  shock  wave 
was  expected  to  be  transformed  first  into  a plastic  wave  and  finally  into 
an  elastic  seismic  wave.  Transient  increases  in  temperature  and  pressure 
were  expected  to  accompany  the  wave  especially  in  fhe  shock  region,  with 
some  permanent  effects  in  all  but  the  elastic  regions. 


The  measurements  initially  proposed  were: 

Part  1:  Measurement  of  pressure  within  the  test  chamber  after  the 
shot.  Since  one  of  the  questions  to  be  answered  by  this  experiment  was 
whether  radioactivity  would  penetrate  to  the  surface  after  the  shot,  and 
since  an  elevated  pressure  in  the  shot  cavity  would  be  the  most  obvious 
means  of  driving  radioactive  gas  or  fine  particulate  material  to  the  sur- 
face, the  measurement  of  pressure  in  the  cavity  was  considered  important. 
It  was  proposed  to  drill  three  very  narrow  vertical  shafts  into  the  region 
near  the  test  chamber  before  the  shot.  While  the  interior  portion  of  the 
holes  might  be  eliminated  by  the  explosion,  it  was  expected  that  these  holes 
would  furnish  a convenient  means  of  inserting  pressure  probes  into  the  cav- 
ity after  the  shot.  It  was  also  hoped  that  radiochemical  sampling  would  be 
possible  through  these  access  holes. 

Part  2:  Measurement  of  temperature  of  the  unaltered  rock.  Eight 
channels  of  temperature  measurement  were  to  be  provided.  Probes  were 
to  be  inserted  in  bare  holes  before  the  shot.  It  was  intended  to  measure 
the  long-term  temperature  changes  due  to  heat  transfer  before,  and  if  pos- 
sible, following  the  explosion. 


Part  3:  Measurement  of  temperature  of  the  molten  rock  region.  A 
similar  system  providing  four  channels  of  temperature  measurement  was 
proposed  for  measuring  the  temperature  of  the  molten  rock  expected  to  be 
contained  in  the  cavity.  It  was  expected  that  the  access  holes  mentioned 
in  Part  1 might  be  used  for  the  insertion  of  temperature  sensing  probes 
after  the  shot. 


Better  knowledge  of  the  expected  results  of  the  shot,  which  became 
available  soon  after  the  starting  date  of  April  24,  1957,  caused  a modifi- 
cation of  the  program. 


The  rocks  composing  the  mesa  in  which  the  explosion  was  to  take 
place  consist  of  rhyolites  and  other  volcanic  tuffs,  with  measurable  air 
space  wvter  content.  A preliminary  shot,  using  50  tons  of  high  ex- 
plosive, showed  that  there  would  be  a considerable  degree  of  fracturing 
and  dislocation  of  the  rock  in  the  vicinity  of  the  explosion  cavity.  The 
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original  concept  of  the  containment  of  the  burst  in  a sealed  cavity,  with 
relative  stability  of  other  regions,  had  to  be  modified.  It  now  seemed  likely 
that  the  cavity  would  be  in  contact  with  fissures,  and  that  any  overpressure 
which  might  originally  be  present  would  very  soon  be  relieved. 

At  the  same  time,  calculations  were  started  which  were  designed  to 
estimate  the  extent  of  heat  conduction  through  the  rock.  Even  the  first 
rough  results  showed  clearly  that  measurable  temperature  increases  due  to 
conduction  would  be  confined  to  distances  of  a few  feet  from  heat  sources 
for  periods  of  months  or  even  years.  It  therefore  appeared  that  probes 
would  have  to  be  installed  in  the  immediate  vicinity  of  the  test  chamber, 
and  (if  installed  before  the  shot)  would  probably  not  survive  the  explosion. 

In  view  of  these  results  the  program  of  investigation  was  changed. 

Part  1 (pressure  measurement)  was  greatly  reduced  in  scope.  Pressure 
measurement  in  the  cavity  was  abandoned,  and  it  was  decided  to  confine 
pressure  measurements  to  gauge  installations  at  the  blast  doors  of  the  tun- 
nel. These  gauges  were  to  indicate  pressures  on  the  inside  to  observers  on 
the  outside  of  the  blast  doors. 

Parts  2 and  3 (temperature  measurement)  were  combined.  All  static 
temperature  tests  were  to  be  made  by  means  of  probes  inserted  in  holes 
which  were  to  be  drilled  after  the  shot.  At  this  stage  it  was  intended  that 
all  holes  would  be  drilled  in  a vertical  direction  from  the  top  of  the  mesa. 

One  hole  was  to  penetrate  the  cavity,  and  temperatures  in  it  were  to  be 
logged  by  a high-temperature  probe.  Other  holes  were  to  be  drilled  to  pre- 
determined points  near  the  cavity,  permanent  probes  were  to  be  installed, 
and  the  change  of  temperature  with  time  was  to  be  measured. 

In  addition  it  was  proposed  that  four  to  six  probes  be  developed  and 
installed  which  would  measure  the  transient  temperature  rise  in  the  rock  due 
to  the  passage  of  the  shock  wave.  These  measurements  were  to  be  made  at 
distances  of  25  to  500  ft  directly  above  the  test  chamber.  It  was  hoped  that 
these  measurements  would  help  to  define  the  behavior  properties  of  the  rock 
in  the  shock  state.  Because  of  the  short  time  scale  it  was  agreed  that  re- 
cording apparatus  furnished  by  the  Sandia  Corporation  would  be  utilised. 

1.3  THEORY 

1.  3. 1 Conduction  of  Heat  and  Temperature  Distribution  in  Tufaceous  Material. 

The  following  calculations  were  carried  out  in  order  to  determine  the 
nature  of  the  temperature  distributions  expected  after  the  shot,  and  to  help 
in  planning  the  measurement  program. 

In  addition,  it  was  hoped  that  these  calculations  might  be  useful  after 
the  shot  in  attempting  to  determine  the  physical  form  and  distribution  of  the 
energy  deposited  by  the  bomb  in  the  tuff.  These  calculations  would  be  use- 
ful also  in  normalising  the  total  heat  energy  present  (giving  an  indication 
of  the  bomb  yield)  by  comparison  with  the  calculations  for  a specified  hea* 
source.  The  total  energy  present  would  be  given  by 

Q = j*4ir  Cd  v(r)  r2  dr 


where  v(r)  is  the  experimentally  determined  temperature  distribution,  C 
is  the  specific  heat  of  the  rc<  k,  d its  density  and  r is  the  radius.  By  com- 
paring v(r)  at  convenient  r with  the  calculations  it  was  hoped  that  the  ex- 
perimental curves  could  give  indications  as  to  the  physical  configuration  of 
hot  material  at  very  small  radii. 

It  should  perhaps  be  emphasised  that  all  results  in  this  section  deal 
only  with  conduction  of  heat  through  a homogeneous  rock.  Experience  in- 
dicates that  additionaT'methods  of  heat  transfer  may  be  operative,  especially 
convection  through  fractures  and  fissures  and  the  transfer  of  heat  by  means 
of  vapors  of  volatile  substances  which  may  recondense  in  cooler  regions. 

If  this  occurs  these  calculations,  as  presented  below,  are  not  representa- 
tive of  the  true  situation. 

The  calculations  were  all  made  from  a single  set  of  physical  assump- 
tions, based  partially  on  the  data  regarding  rock  characteristics  which  were 
available  at  the  time  the  calculations  were  performed.  Subsequently  better 
data  have  become  available;  however,  these  data  do  not  alter  the  qualitative 
aspects  of  the  results  and  it  was  not  felt  worthwhile  to  repeat  the  calcula- 
tions using  the  newer  data.  The  following  assumptions  were  made: 

(a)  Heat  energy  liberated,  Q,  is  1.8  X 10^  cal. 

(b)  Rock  contains  9%  HoQ  and  the  heat  required  to  raise  it  to  2000*C 
is  728  cal/g. 

(c)  Density  of  rock,  d,  is  taken  as  1.8Sf  g/cm3.  Air  space  taken 
as  16%. 

, (d)  Heat  conductivity  of  rock,  K,  is  taken  as  1.  7 X 10” 3 cal  cm-1 
•C  sec-1. 

(e)  Heat  capacity,  C,  taken  as  0.  330  cal/g  *C. 

(f)  This  gives  for  the  diffusivity  S = K/Cd  the  value  2.66  X 10”  3. 

For  comparison  we  cite  the  following  corresponding  values  for  four  other 
types  of  earth  material: 


Material 

_d 

C 

K 

S 

"Average  rock" 

4.  2X10” 3 

11.8X10“ 

Brick 

2.6 

0.20 

2.  0X10”3 

3.  8X10" 

Soil  (average) 

2.5 

0.2 

2.  3X10" 3 

4.6X10” 

Soil  (sandy,  8%  moist) 

1.75 

0.24 

1.4X10” 3 

3.  3X10” 

(It  will  be  seen  that  the  tufaceous  rock  of  the  mesa  is  presumed  to 
be  closer  to  brick  or  sandy  soil  than  to  average  rock,  because  of  its  rela- 
tively open  structure.) 

(g)  In  some  calculations  it  is  assumed  that  a volume  of  rock  has  been 
melted  and  heated  to  2000*C  by  the  explosion,  and  that  it  fills  a spherical 
cavity.  Calculations  based  on  the  data  given  above  yield  678  cm  for  a,  the 
radius  of  such  a cavity. 

In  all  succeeding  calculations  quantity  v will  designate  the  temperature 
difference,  (T-Tfl)  (v=0  before  shot  throughout  mesa).  For  a homogeneous 
isotropic  solid,  such  as  we  assume  the  mesa  to  be,  the  differential  equation 
of  heat  conduction  is 


assuming  that  the  diffusivity  is  constant  over  the  time  and  temperature 
ranges  under  consideration*  and  we  can  investigate  solutions  corresponding 
to  different  boundary  conditions. 


1.  3. 1. 1 Point  Source.  The  simplest  case  to  consider  is  that  in  which 
the  entire  quantity  of  heat  is  liberated  instantaneously  at  a single  p--nt.  In 
this  case  the  equation  is  satisfied  by: 

Q -(rZ/4St) 

8Cd(sSt)i/r2 


V = 


Putting  in  the  assumed  figures*  this  becomes 

v = 4.  72X1014  t"3/2  exp  {-r2/0.  01064t) 

where  time  t is  in  seconds  and  radius  r in  centimeters.  The  values  for  v 
obtained  in  this  manner  are  given  in  Table  1. 1. 


It  can  be  seen  from  the  first  line  of  Table  1.  1 that  at  any  given  point 
the  temperature  rises  to  a maximum  and  then  drops.  The  maximum  tem- 
perature and  the  time  of  its  occurrence  can  be  calculated  as  follows: 

v = A t"3/2  exp  (-B/t); 


A 

dv 

St 


4. 72X10 14;  B = r 2/0. 01064 
A e'®/*  (Bt-7/2  - -If5/2)  = 0 

62.  7 r2;  v = 2. 125X101 1 f3 
max 


Values  of  t and  v for  given  r are  given  in  Table  1.  2. 
max  max  ^ 

Keeping  in  mind  that  one  year  is  3.  16x10  seconds,  we  see  that  con- 
siderable times  will  elapse  before  any  temperature  rise  can  be  detected  at 
radial  distances  for  which  this  approximation  might  be  valid. 

1.  3. 1.  2 Homogeneous  Spherical  Source.  In  order  to  improve  the  calcu- 
lation, we  examine  some  other  solutions  which  may  give  a better  approxi- 
mation to  reality.  One  such  calculation  lets  the  heat  be  evolved  instantan- 
eously and  homogeneously  inside  a spherical  volume  of  radius  a. 

Carslaw  and  Jaeger  give  a partial  solution  for  this  case  (Reference  1): 


v s 


2r(irSt)1/Z 


e-(r-r')2/4St 


dr* 


where  V is  the  initial  temperature  inside  the  sphere.  This  can  be  integrated, 
and  making  the  substitution  V = 3Q/4ea^Cd,  we  obtain: 

. — a, — . .w/®'! . wr «rf 

8 r a3  Cd  \ L J \ 2V51 

The  solution  of  the  integral  is  recorded  in  Appendix  A. 

In  the  calculation  the  same  constants  as  before  are  used,  with  the  addi- 
tion of  a a 678  cm.  The  results  are  given  in  Table  1.  3. 

Since  a = 678  cm,  values  for  r = 500  cm  are  omitted  as  falling  inside 
the  shell,  where  latent  heats  of  fusion  complicate  the  picture.  It  will  be 
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TABLE  1. 1 TEMPERATURE  RISE  v IN  *C  - INSTANTANEOUS 


POINT  SOURCE 


V*'  *ec  6 

\ 1 X 10° 

r,  cm  N 

5 X 106 

1 X 107 

2 X 107 

3 X 107 

1 X 108 

500 

2.95X10"5 

38.4 

1424  1630 

1313 

373 

1000 

7.  17X10" 36 

2.90X1  O’4 

1.24 

48.0 

125 

184 

1500 

6.  84X10'®7 

1.81X10”14 

9.78X10"6 

.135 

2.49 

57  X) 

2000 

< io- : 100 

9.37X10"29 

7.04X10"13 

3.62X10”5 

0.0104 

110 

2500 

< 10" 100 

4.02X10"47 

4.71X10"22 

9.27X10"10 

9.01x1 0"4 

1.33 

3000 

<10" 100 

1.43X10"69 

2.75X10"33 

2.26X10"15 

1.63x10”* 

0J00 

TABLE  1.2  v and  t - INSTANTANEOUS  POINT  SOURCE 
max  max 


500 

1.57X107 

1700 

600 

2.  26x  » 

984 

700 

3.08X  •' 

618 

800 

4.02X  » 

414 

900 

5.08X  " 

291 

1000 

6.27X  »» 

212 

1200 

9.03X  " 
1.23X10° 

123 

1400 

77.4 

1600 

1.61X  •» 

51.9 

1800 

2.03X  •• 

36.4 

2000 

2.51X  »' 

26.6 

2200 

3.03X  " 

20.0 

2400 

3.61X  »• 

15.4 

2600 

4.  23X  '• 

12.  1 

2800 

4.92X  »' 

9.68 

r,  cm 

t i sec 

max 

vmax* 

3000 

5. 64X1 08 

7.  87 

3300 

6.83X  •» 

5.91 

3600 

8. 13X  »' 

4.55  ' 

3400 

9.54X  •• 

3. 58 

4200 

1.11X109 

2.87 

4500 

1.27X  •• 

2.  33 

4800 

1.44X  »» 

1.92 

5100 

1.63X  •' 

1.60 

5400 

1.83X  » 

1.35 

5700 

2.04X  •• 

1.15 

6000 

2.  26X  " 

0.984 

6500 

2.65X  « 

0.774 

7000 

3.08X  »' 

0.618 

TABLE  1.  3 TEMPERATURE  RISE  v IN  *C  - HOMOGENEOUS 


SPHERE  SOURCE 


\ t,  sec 

r,  cms 

1X106 

5X106 

1X107 

2X107 

3X107 

1X108 

1000 

0.008 

33.0 

103 

185 

210 

167 

1500 

< 10“  3 

< 10“3 

< 10"3 

4.94 

15.5 

61.7 

2000 

< 10“3 

< 10"3 

< 10"3 

< 10“  3 

0.  300 

15.2 

2500 

< 10“3 

< 10"3 

< 10“  3 

< 10“3 

0.007 

2.  70 

3000 

< 10“3 

< 10“  3 

< 10"3 

< 10“3 

< 10“3 

0.253 
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seen  that  for  large  values  of  r and  t the  results  tend  to  approach  those  in 
the  previous  table.  Because  of  limitations  in  the  available  tabulations  of 
the  erf  function,  no  values  of  v under  0.  001  *C  are  given. 

1.  3.  1.  3 Spherical  Cavity  Surface  of  Constant  Temperature.  Another 
possible  approximation. is  given  by  assuming  that  heat  proceeds  from  the 
surface  of  a sphere  of  radius  a into  an  otherwise  unbounded  homogeneous 
region,  the  surface  of  the  sphere  being  kept  at  a constant  temperature  V. 
For  this,  the  temperature  rise  at  any  time  t and  radius  r (>  a)  is  given  by 


v = 


aV 

r 


erfc 


The  use  of  this  approximation  implies  a supply  of  heat  to  the  surface  of 
the  cavity,  e.  g. , by  freezing  of  molten  rock. 

We  use  the  same  constants  as  before,  letting  V = 2000*C,  and  we 
obtain  results  which  are  given  in  Table  1.  4. 

As  might  be  expected,  this  approximation  gives  no  maximum  v at 
any  r.  Because  of  the  implied  continuous  generation  of  heat,  this  equation 
would  not  be  expected  to  be  valid  at  large  values  of  r and  t.  However,  it 
might  be  the  best  approximation  for  small  values  of  t and  (r-a),  if  our  phys- 
ical assumptions  are  valid  or  if  better  ones  can  be  found. 

TABLE  1.  4 TEMPERATURE  RISE  v IN  *C  - SPHERICAL  SURFACE 
AT  CONSTANT  TEMPERATURE 


v t,  sec 

r,  crn^ 

1 x io6 

5 X 106 

1 x io7 

2 X 107 

3 X 107 

1 X 108 

1000 

0.014 

65.6 

221 

439 

570 

894 

1500 

< 10"3 

< 10*3 

< IO-3 

10.6 

35.8 

235 

2000 

< 10’3 

< IO-3 

< 10"3 

< 10"3 

0.  637 

47.4 

2500 

< 10*3 

< 10"3 

< 10‘3 

< 10’3 

0.005 

6.  89 

3000 

< 10"3 

< 10"3 

< 10"3 

< 10° 

< 10° 

0.655 

I.  3.  1.  4 Heat  Generation  at  the  Spherical  Surface.  Yet  another  approx- 
imation can  be  made  if  we  assume  that  the  heat  is  generated  instantaneously 
at  the  surface  of  a spherical  shell  of  radius  a,  surrounded  inside  and  outside 
by  homogeneous  rock.  For  this  approximation  we  have  the  equation 


8irraCd(wSt) 

Making  the  same  assumptions  as  before,  we  obtain  temperature  rises 
from  this  equation  which  are  shown  in  Table  1.  5.  Since  the  error  function 
does  not  occur  here,  small  values  of  v can  be  more  easily  tabulated. 
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TABLE  1.5  TEMPERATURE  RISE  v IN  *C  - SPHERICAL  SURFACE  SOURCE 


r,  cm 

• *«c  A 

1 X 10° 

5 X 106 

1 X 107  2 X 107 

3 X 107 

1 X 108 

1000 

0.109  118 

221  249  244 

155 

1500 

3.  3X10" 25 

1.68X10"3 

0.682  18.5 

27.2 

64.0 

2000 

4.  3X10"69 

2.  2X1 O' 12 

2. 15X10'5  0.0562 

0.  709 

17.8 

2500 

< 10"100 

2. 6X10“ 25 

6. 6X10" 12  2.  78X10"5 

4. 12X10" 3 

3.27 

3000 

< io-100 

2.  7X10"42 

1.9X10"20  1.4X10"9 

5.2X10"6 

0.389 

1.3. 1.5  Instantaneous  Source  at  the  Surface  of  a Spherical  Cavity.  A 
final  approximation  is  the  one  given  in  Reference  1,  page  303.  Here  we 
again  have  an  instaneous  spherical  surface  source  at  r • a,  but  we  also  as- 
sume aero  heat  flux  across  the  boundary  r = a.  This  might  be  considered 
equivalent  to  an  empty  cavity  whose  walls  are  heated  by  radiation.  For  this 
we  write  the  formula: 


• *£*{(ir),/2  - ! - f?  * fi]  - ( 


r-a 

ZsfSt1 


Using  the  same  assumptions  as  before,  this  yields  the  results  for  v 
given  in  Table  1.6.  Plots  of  results  obtained  from  these  and  the  previous 
approximations  are  given  in  Fige.  1. 1,  1.  2 and  1.  3,  and  show  a general 
correspondence  in  trend. 


TABLE  1. 6 TEMPERATURE  RISE  v IN  *C  - SPHERICAL  CAVITY 


SURFACE  SOURCE 


r.  cm' 

sec  c 

1 X 10° 

5 X 106 

1 X 107 

2 X 107 

3 X 107 

1 X 10® 

1000 

0.211 

206 

363 

360 

325 

147 

1500 

< 10"3 

< 10"  3 

1.23 

19.2 

42.5 

73.9 

2000 

< 10"3 

< 10"3 

< 10"3 

0.099 

1.20 

23.2 

•> 

2500 

< 10"3 

< io"3 

< 10"  3 

< 10"3 

< 10"3 

4.55 

3000 

< 10"3 

< 10"3 

< 10"3 

< 10"3 

< 10"3 

0.58 
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1.  32  TEMPERATURE  RISE  IN  TUFF  DURING  PASSAGE  OF  SHOCK  WAVE 

Calculations  of  the  expected  temperature  rise  in  tuff  may  be  made  in 
those  regions  in  which  a shock  wave  will  be  supported.  These  calculations 
as  reported  here  follow  the  work  of  previous  investigators  (Reference  3). 

We  assume  that  the  tuff  can  be  represented  by  an  idealised  mixture  of 
incompressible  solid  and  air,  and  an  ideal  gas.  Let  4 = P8/P0  = pressure 
ratio  at  the  shock  front,  T = temperature,  p = density,  and  P = pressure. 
Subscripts  o,  f,  and  s refer  to  initial,  final,  and  shock  conditions,  respec- 
tively. 

The  Rankine-Hugoniot  adiabat  yields,  for  an  ideal  gas 


£ 


Let  mg  be  the  mass  per  unit  volume  of  air  and  mj  be  the  mass  per  unit 
volume  of  the  solid.  Then  the  increase  in  heat  content  of  the  air  is 

Q.lr  = “2  VT'T.>  ' ”2  • 

where  h^  is  the  specific  heat  of  the  air. 

Assuming  very  rapid  mixing  of  the  air  and  incompressible  solid  so 
that  equilibrium  heat  transfer  occurs,  the  total  heat  content  increase  for 
the  mixture  just  equals  that  for  the  air  so  that 


or 


whence 


' Qmi, 


(m,  1>,  + m2  h2)  (T,  - T0>  * i»2  1>2  Ifejl 

£ 


* i + 


m2  h2 


mi  hi  + mg  h2 


Since  for  most  cases  h,  “ h_  , 

Calculation  of  the  temperature  rise  at  any  radius  thus  resolves  itself 
to  a determination  of  |,  the  pressure  ratio  as  a function  of  radius. 

Several  methods  are  available  for  estimating  the  peak  shock  pressure 
at  various  radii  from  ground  sero.  The  most  straight-forward  of  these 
methods  involves  the  direct  machine  integration  of  the  hydrodynamic  conser- 
vation equations  forward  in  time  for  each  radius.  The  difficulties  inherent 
in  such  calculations  are  great  and  the  question  of  initial  conditions  involving 
wave  shapes  is  troublesome.  Nevertheless,  calculations  have  been  made 
in  this  fashion  for  some  media. 


In  addition,  several  approximate  schemes  have  been  developed  for 
dealing  with  the  propagation  of  shock  waves  in  media.  Some  of  these  are 
the  Penny-Das  Gupta,  Kirkwood-Bethe,  Kirkwood-Brinkley  (referred  to 
in  Reference  2),  Porsel  (Reference  3),  etc.  schemes.  One  derived  curve 
of  pressure  versus  distance  derived  by  UCRL  is  given  in  Fig.  1.4.  Using 
this  data  the  expected  temperature  rise  as  a function  of  distance  is  plotted 
in  Fig.  1.5.  The  data  used  in  calculating  these  figures  are: 
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CHAPTER  2 
PROCEDURE 


The  construction,  calibration  and  field  installation  of  equipment  used 
to  implement  the  three  parts  described  in  Section  1.  2 are  discussed  in  this 
chapter. 

2. 1 BLAST  DOOR  INSTRUMENTATION 

The  access  tunnels  to  the  site  of  the  nuclear  shot  and  to  the  sites  of 
preliminary  shots  were  closed  by  heavy  blast  doors  set  in  concrete  frames. 

It  was  desirable  to  ascertain  after  the  shot  the  pressure  and  composition  of 
the  air  or  gas  behind  these  blast  doors,  and  Broadview  Research  Corporation 
undertook  the  task  of  providing  instrumentation  for  this  purpose. 

The  technique  used  here  consisted  of  pipes  running  through  the  con- 
crete frames  of  the  blast  doors  and  connected  to  the  gage  stations  with  cop- 
per tubing.  The  assembly  started  with  one-inch  iron  pipe  sections,  three 
feet  long,  which  were  inserted  in  holes  drilled  through  the  blast  doors  and 
then  grouted  in.  These  were  connected  by  standard  pipe  couplings  to  3/8- 
inch  (nominal)  type-K  soft-copper  water  tubing  (0.  500-in.  o. d. , 0.049-in. 
wall  thickness).  This  ran  to  the  gage  stations,  passing  through  any  additional 
blast  doors  by  means  of  one-quarter-inch  iron  pipe,  threaded  at  both  ends 
and  coupled  to  the  copper  tubing.  Three  kinds  of  gage  stations  were  used. 

The  first  type  (No.  1)  consisted  of  a copper  tubing  cross  inserted  into  the  cop- 
per tubing  assembly. 

Two  openings  contained,  respectively,  (a)  a needle  valve  to  the  tunnel 
and  (b)  a needle  valve  and  pressure  gage.  One  of  these  stations  was  in- 
serted in  each  line  on  the  outside  of  each  blast  door.  The  second  type  (No.  2) 
was  installed  at  the  tunnel  exits  and  was  identical  with  No.  1 except  that  a 
T was  used  instead  of  a cross  at  the  end  of  the  line.  The  third  type  (No.  3) 
was  a gas-and-dust  sampling  station.  One  of  these  stations  was  installed  at 
the  end  of  the  one-inch  pipe  through  the  originating  door  at  each  assembly. 

It  consisted  of  a gate  valve  and  sampling  head.  Gate  valves  were  chosen  in 
order  to  provide  a straight-line  path  to  the  sample  head  and  to  minimize  settling 
out  of  dust. 

Pressure  gages  were  by  Ashcroft,  type  1379a.  They  had  ranges  of 
0-160  psig  and  were  calibrated  to  0.5%  accuracy. 

Three  tunnels  were  instrumented  in  this  manner  by  Broadview  per- 
sonnel. Tunnels  U-12-B  and  U-12-D  were  instrumented  in  July,  1957,  and 
tunnel  U-12-C,  in  September,  1957.  Two  blast  doors  were  instrumented  in 
U-12-B,  two  in  U-12-C  and  one  in  U-12-D.  Approximately  2,  700  ft  of  cop- 
per tubing,  twelve  pressure  gages  and  the  necessary  auxiliary  equipment 
(pipe,  couplings,  etc.)  were  used.  The  outer  ends  of  the  assemblies,  those 
with  the  No.  2 terminators,  were  installed  in  rugged  crypts  near  the  tunnel 
entrances.  After  the  systems  were  assembled,  they  were  pressure-tested  at 
1000  psi  with  the  gage  valves  closed,  and  all  leaks  were  eliminated.  A final 
test,  including  the  gages,  was  run  at  150  psi,  and  the  equipment  was  left  in 
place  with  the  gage  valves  shut. 

2.  2 INSTRUMENTATION  FOR  MEASUREMENT  OF  TRANSIENT  TEMPERA- 
TURE RISE  ACROSS  SHOCK  WAVE  IN  TUFF 

The  requirements  which  must  be  fulfilled  by  a sensor  capable  of  re- 
producing the  transient  temperature  rise  across  a shock  wave  are: 


| 
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1.  Very  rapid  response  (10  to  10  seconds). 

2.  Electrical  output  signal  for  remote  recording. 

3.  Intimate  coupling  to  sampled  medium. 

4.  Insensitivity  to  stimuli  other  than  temperature. 

Initially,  it  might  seem  that  either  a thermocouple  or  a resistance 
thermometer  might  be  appropriate  for  the  designated  purpose.  However, 
it  was  soon  seen  that  conventional  thermocouples  could  not  give  a sufficiently 
rapid  response.  In  conventional  thermocouples  there  is  a welded  bead  or 
other  high-mass  irregularity  at  the  metal  junction.  Thin-sheet  thermo- 
couples have  been  plated  onto  metallic  bases,  but  this  method  was  not  appli- 
cable in  this  case.  It  was  therefore  decided  to  develop  a resistance  thermom- 
eter which  had  to  be  small  in  order  to  allow  rapid  heating.  It  also  had  to 
be  grouted  into  the  rock  in  order  to  couple  with  the  rock  well  enough  to  make 
sure  that  the  temperature  rise  measured  would  be  that  of  the  rock  and  not  of 
the  thermometer.  This  immediately  raises  a problem.  Without  a grouting 
or  binding  material,  contact  with  the  rock  cannot  be  adequate;  but  if  such  a 
material  interposes  between  the  thermometer  and  the  rock,  its  own  tempera- 
ture will  be  the  one  recorded.  It  was  thus  necessary  to  choose  a grouting 
material  which  was  similar  in  its  properties  to  the  rock.  The  material  se- 
lected was  a product  of  U.S.  Gypsum  Company,  Hydro-Stone,  which  is  simi- 
lar to  plaster  of  Paris  in  appearance  and  properties.  The  density  and  fraction 
of  air  space  of  Hydro-Stone  aggregates  averaged  around  1.  70  and  0.  25,  re- 
spectively. All  measurements  of  transient  temperature  therefore  relate  to 
an  environment  of  Hydro-Stone. 

The  resistance  of  metals  also  changes  as  a result  of  strain,  and  it 
was  important  to  choose  a thermometer  design  for  which  this  effect  would 
be  small.  This  implies  a high  change  in  resistance  with  temperature.  Suit- 
able metals  are  platinum,  copper,  and  gold  having  temperature  coefficients 
of  resistivity  of  3-4  X 10"  3 *C_1. 

Three  approaches  to  this  problem  were  tried.  They  were: 

1.  "Platinum-bright"  strips. 

2.  Platinum  or  copper  wire. 

3.  Gold  foil. 

Liquid  bright  platinum  is  a product  of  Hanovia  Company  of  Newark, 

N.J.,  and  consists  of  a suspension  or  solution  of  platinum  salts  in  an  organic 
vehicle.  Upon  heating,  this  material  is  reduced  to  platinum.  When  applied 
1:  '•lass  or  other  smooth  surfaces,  platinum  is  obtained  as  a smooth  coherent 
layer.  Such  layers  have  been  successfully  used  as  transient  resistance  ther- 
mometers by  Rabinowics,  Jessey,  and  Bartsch  of  the  Guggenheim  Aeronaut- 
ical Laboratory  of  the  California  Institute  of  Technology,  and  by  Dr.  D. 
Bershader  and  his  group  at  Lockheed  Corporation.  The  attempt,  here,  was 
made  to  form  suitable  conducting  layers  of  platinum  on  Hydro-Stone  surfaces. 
Smooth  Hydro-Stone  surfaces  were  obtained  by  casting  on  glass.  Lines  were 


then  drawn  on  small  (0. 8-1  cm)  shards  of  this  material  and  the  platinum  was  l 

reduced  by  ignition  at  600*C.  Pieces  which  gave  suitable  resistances  (in  the 

vicinity  of  300  ohm)  were  selected  and  small  drops  of  solder  were  applied  to 

each  end  of  the  platinum  line.  With  suitable  care  these  could  be  made  to 

wet  the  platinum  surface.  Copper  wires  were  then  set  into  the  solder  beads,  » 

and  the  entire  assembly  was  put  into  a mold  which  was  filled  with  Hydro- 

Stone.  About  eighty  of  these  elements  were  prepared.  Water -proofing  with 

paraffin  and  with  collodion  proved  insufficiently  effective,  but  two  coats  of 

Epon  820  epoxy  resin,  mixed  with  Epon  curing  agent  U in  4: 1 ratio,  gave 

- 20- 


adequate  protection  from  moisture.  It  was  estimated  that  the  thickness  of 
the  conducting  platinum  layer  was  0.  5-1  p.  Since  this  method  gave  excellent 
bonding  to  the  Hydro-Stone  and  negligible  thermometer  lag,  it  was  hoped  that 
it  could  be  adopted.  Unfortunately,  thermometers  which  were  prepared  by 
this  method  showed  fluctuations  in  resistance;  in  some  cases  all  conductivity 
was  lost.  The  platinum  layer  had  no  tensile  strength  except  that  of  the  Hydro- 
Stone  substrata  and  it  appeared  that  stresses  in  the  Hydro-Stone  cause 
cracks  (even  long  after  casting)  which  could  partially  or  totally  sever  the  con- 
tinuity of  the  platinum  strip.  This  method  was,  therefore,  eliminated  from 
consideration. 

Wire  resistors  were  produced  from  0.001-in.  platinum  wire,  obtained 
from  Wildberg  Brothers  Smelting  and  Refinery  Company  of  San  Francisco, 
and  from  0.002-in.  General  Electric  copper  magnet  wire,  covered  with 
Formex  insulation.  The  platinum  wire  resistors  required  25  ft  of  wire  and 
were  wound  on  cylindrical  l/2  in.  x 1/2  in.  Hydro-Stone  cores,  especially 
grooved  to  keep  the  wires  from  touching.  Small  dabs  of  Lix-Crax  ceramic 
cement  were  used  to  keep  the  wire  from  unwinding.  The  copper  resistors 
were  wound  on  1 in.  x 1-1/2  in.  cards.  One -hundred  feet  of  wire  were  re- 
quired. The  resistors  in  each  case  were  set  in  molds  and  buried  in  cast 
Hydro-Stone.  Copper  lead  wires  were  soldered  on  (in  the  case  of  the  copper 
wire). 

The  small-scale  tests  (described  later)  which  were  run  on  these 
thermometers  showed  large  unexplained  signals,  which  were  attributed  to 
induction  effects.  Non-inductively  wound  copper  resistors  showed  some  im- 
provement. In  the  end,  however,  it  was  decided  not  to  utilize  these  resistors 
because  of  their  significantly  greater  mass  and  poorer  bonding  to  the  Hydro- 
Stone. 

The  gold  leaf  resistors  which  were  finally  selected  were  made  from 
gold  leaf  manufactured  by  Witt  Coe  Manufacturing  Company  of  Providence, 

R.I.  This  leaf  came  in  sheets  3-3/8  in.  square.  Sheets  of  gold  leaf  were 
placed  on  the  wetted  surface  of  Hydro-Stone  slabs,  1/8  inr  thick,  which  had 
been  previously  cast  against  clean  glass.  The  sheets  of  gold  leaf  were  picked 
up  by  one  edge,  using  a moistened  enameled  pencil  to  which  the  gold  adhered 
to  form  a small  flag.  The  free  end  of  the  leaf  was  brought  in  contact  with  the 
wet  Hydro-Stone.  Interfacial  tension  between  the  gold  and  the  water  then  drew 
the  gold  leaf  down  to  the  wet  surface  of  the  Hydro-Stone.  This  procedure, 
performed  in  draft-free  room,  produced  smooth  gold  films  on  the  Hydro- 
Stone  after  a small  amount  of  practice.  After  the  Hydro-Stone  dried,  the  gold 
leaf  adhered  to  it.  These  sheets  of  gold-coated  Hydro-Stone  were  sawed 
into  pieces  l/2  in.  square.  The  surfaces  of  the  squares  were  scribed  to 
produce  the  resistance  elements.  A piece  of  hypodermic  tubing,  0.025-in. 
o.  d. , was  filed  to  produce  a square  end  beveled  on  the  inside  to  produce  a 
sharpened  edge.  The  tube  was  then  bent  into  a V,  and  the  unsharpened  end 
held  in  a pin  vise.  When  this  small  tubular  chisel  was  drawn  over  the  sur- 
face of  gold  leaf  on  Hydro-Stone,  it  cut  a smooth  groove  in  the  surface  with- 
out appreciable  tearing  of  the  adjacent  gold  leaf.  Equally  spaced  grooves 
were  cut  in  the  gold-filmed  surface  of  the  Hydro-Stone  in  such  a pattern  that 
an  extended  electrical  path  was  produced  in  the  gold  leaf.  The  spacing  of  the 
grooves  was  adjusted  to  produce  a resistance  of  approximately  300  ohms 
across  the  element.  The  0.010-in.  copper  leads  were  cemented  to  the  back 
of  the  element  in  such  a way  that  their  ends  were  adjacent  to  opposite  ends  of 
the  gold  leaf  element.  Connection  to  the  gold  leaf  was  made  by  conducting 
silver  paint.  The  element  leads  were  soldered  to  connecting  cables  and  the 
assembly  centered  in  a greased  aluminum  mold  which  was  subsequently  filled 
with  Hydro-Stone  to  produce  the  probe.  The  probes  were  waterproofed  with 
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two  coats  of  Epon  820  resin  as  above  and  were  then  copper  plated  to  a thick- 
ness of  approximately  0.005  in.  This  copper  plating  was  designed  to  com- 
plete electrical  shielding  of  the  probe  and  was  connected  to  the  cable  shield, 
which  was  used  as  a current  lead.  For  further  waterproofing,  a plast:  seal 
was  cast  over  the  cable  end  of  the  probe  and  the  first  1-1/2  in.  of  the  t^jle. 

The  plastic  used  was  Scotchcast  No.  4,  and  a 3M  liquid  rubber-type  sealant 
was  used  to  bond  the  Scotchcast  to  the  plastic  cable  covering.  Molds  were 
L machined  from  Teflon  bar  stock;  fairly  large  pouring  holes  were  found  neces- 

sary. Figures  2.  1 and  2.  2 show  these  probes  in  various  stages  of  production. 

In  order  to  test  the  response  of  the  transient  temperature  probes,  . 
scaled  explosive  tests  were  undertaken.  The  scaling  was  on  the  basis  W*'  3 d * 
const,  where  W is  the  TNT  equivalent  of  the  charge  and  d is  the  distance  from 
the  center  of  the  charge  to  the  sensing  element  of  the  probe.  The  charges 
used  were  Composition  A 110-gram  spheres.  Since  these  charges  are  about 
2 in.  in  diameter,  and  d values  ranged  from  6 to  24  in. , it  is  evident  that 
great  precision  is  impossible;  fortunately  it  was  not  necessary  for  the  quali- 
tative results  desired.  Experiments  were  carried  out  in  Hydro-Stone  and  oil 
media.  For  the  Hydro-Stone  experiments,  the  probes  were  embedded  in  cy- 
lindrical Hydro-Stone  blocks  9 to  15  in.  in  diameter  depending  on  d.  The 
probes  were  buried  in  at  least  6 in.  of  Hydro-Stone,  with  leads  coming  through 
one  flat  surface  of  the  cylinder,  and  the  charge  was  set  in  a hemispherical 
cavity  on  the  other  flat  surface.  For  the  oil  tests,  the  probes  were  put  in  a 
cavity  in  the  Hydro-Stone  block  which  was  filled  with  lubricating  oil.  Charge 
and  block  were  buried  in  sand.  In  most  shots,  control  probes  of  nichrome 
wire  were  placed  next  to  the  test  probe.  Signals  were  obtained  as  oscillo- 
scope traces  and  photographed- 

Twenty-four  shots  were  made.  Results  in  general  showed  sharp  pulses 
upon  arrival  of  the  stress  wave,  with  rapid  trailing  off.  There  were,  however, 
a number  of  unexplained  spurious  peaks,  oscillations  and  signal  failures.  In 
some  instances  these  could  be  explained  by  equipment  failures,  in  others  by 
mechanical  failure  of  the  probe,  inductance  or  strain  effects.  The  gold  leaf 
probes  showed  fewer  unexplained  phenomena  than  others  and  were  selected 
on  this  basis. 

It  should  be  emphasized  that  these  tests  were  very  imperfect  repro- 
ductions of  the  actual  shot  conditions.  In  particular,  the  sealing  laws  do  not 
hold  at  such  close  distances  to  the  charge,  nor  do  they  hold  in  media  other 
than  air.  It  is  evident  from  Section  1.  3.  2 that  in  order  to  obtain  readable 
temperature  rises  we  must  be  in  the  pressure  range  above  10^  bars,  or  in 
the  case  of  conventional  explosives  very  close  to  the  charges  indeed. 

After  the  gold-leaf  probe  design  was  adopted,  about  forty  of  these 
probes  were  manufactured  and  calibrated  using  fixed  temperature  points. 
Twelve  were  sent  to  the  Nevada  Test  Site  for  installation  and  six  were  installed 
at  distances  of  25  to  400  ft  from  ground  zero. 

Records  of  the  signals  from  the  probes  were  taken  on  Sandia  Corpora- 
tion equipment,  six  recording  channels  having  been  made  available  for  these 
measurements  through  their  courtesy.  The  recording  equipment  consisted  of 
a Consolidated  Engineering  Company  System-D  bridge -and- amplifier  system 
whose  output  was  recorded  on  an  Ampex  tape  recorder.  The  overall  frequency 
response  of  this  system  is  less  than  a thousand  cycles  per  second,  much  too 
narrow  for  the  desired  temperature  information,  but  the  best  which  was  avail- 
able within  the  operational  time  scale. 
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2.  3 INSTRUMENTATION  FOR  LONG-TERM  THERMOMETRY 
2.  3. 1«  General 

The  choice  of  an  instrument  for  measuring  long-term  temperature 
changes  in  the  mesa  after  the  shot  lies  at  first  glance  between  thermocouples 
and  resistance  thermometers.  The  thermocouple  has  several  impressive 
advantages.  Thermocouples  are  small,  stable,  wide-range  instruments, 
with  rapid  response  signals  that  are  largely  independent  of  the  length  and 
condition  of  the  connections  between  the  instruments  and  the  recorder,  and 
which  can  be  measured  on  relatively  simple  instruments . However,  they 
are  essentially  low -sensitivity  instruments,  and  it  was  not  anticipated  that 
a precision  of  more  than  0.  1*C  could  be  attained.  Other  less  important  con- 
siderations in  using  thermocouples  were  the  necessity  of  maintaining  a cold 
junction  under  field  conditions  and  the  difficulty  of  maneuvering  with  long 
lengths  of  thermocouple  wire  in  narrow  access  holes. 

Resistance  thermometers  are  based  on  the  change  of  resistance  of  a 
conductor  with  temperature.  Metals  have  small  positive  temperature  coef- 
ficients of  resistance  while  nonmetals  (specifically  semiconductors)  have 
large  negative  temperature  coefficients  of  resistance.  Thus  the  semiconduc- 
tor resistance  thermometers  (thermistors)  are  capable  of  giving  very  great 
precision.  Unfortunately  the  accuracy  of  these  instruments  is  not  equally 
dependable,  as  a certain  amount  of  drift  in  resistance  values  can  sometimes 
be  observed.  It  was  found  that  no  thermistors  were  commercially  available 
which  could  be  depended  upon  to  give  reproducible  results  over  many  months 
under  field  conditions.  A metal  resistance  thermometer  was  therefore  chosen 
for  long-term  temperature  measurements.  Commercially  available  probes 
were  made  of  copper  or  platinum.  Since  it  was  anticipated  that  some  probes 
would  be  left  in  place  for  long  periods  of  time,  stability  was  a very  important 
consideration.  Platinum  resistance  thermometers  were  therefore  selected. 
Platinum  wire  for  thermometer  use  should  pass  two  tests:  (1)  the  ratio  of 
the  resistance  at  100*C  to  that  of  0*C  should  be  1.  392  or  greater,  and  (2) 
the  wire  should  melt  in  a hot  flame  without  sputtering,  thus  indicating  homo- 
geneity and  absence  of  air  inclusions.  Platinum  wire  which  passes  these 
tests  is  used  as  the  primary  standard  of  temperature  over  a range  from 
-183*C  to  +630*C,  which  included  the  expected  values  of  most  of  the  tempera- 
tures to  be  measured,  and  platinum  resistance  thermometers  had  been  suc- 
cessfully used  up  to  1650*C. 

The  planned  measurements  for  this  shot  were  as  follows: 

1.  Measurement  of  rock  temperatures  at  fairly  large  distances 
(lateral)  from  the  site  of  the  explosion.  Access  was  to  be  through  holes 
drilled  from  top  of  mesa.  Probes  were  to  be  left  in  the  holes  for  several 
months . 

2.  Logging  of  vertical  hole  which  was  to  penetrate  into  cavity.  Many 
temperatures  to  be  read  in  relatively  short  time. 

This  plan  was  made  in  the  expectation  that  the  explosion  heat  would 
largely  be  contained  in  a cavity  with  impervious  walls,  and  that  the  main  heat 
transfer  mechanism  would  be  conduction.  In  other  words,  the  calculations 
outlined  in  Chapter  1 were  expected  to  be  fairly  close  to  the  actual  case.  It 
was  therefore  expected  that  temperature  changes  at  distances  of  50-100  ft 
from  gound  zero  would  be  quite  small,  and  that  temperature  probes  for  the 
holes  outside  the  cavity  would  need  to  have  precision  and  accuracy  of  10~2  - 
10"3*C.  On  the  other  hand,  temperatures  in  the  cavity  might  be  over  1000*C, 
and  a rugged  high  temperature  probe  would  be  required  for  logging. 
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The  probes  selected  were  platinum  resistance  thermometers.  The 
resistance  was  100  ohms  at  0*C.  These  were  Weiller-type  instruments, 
encased  in  stainless  steel,  produced  by  Glass-quar  Company  of  New  York 
and  furnished  by  Hallikainen  Instruments  of  Berkeley,  California.  Figure 
2.  3 shows  the  detailed  construction  of  these  probes. 

Twelve  probes  were  originally  ordered  for  measurements  at  tempera- 
tures up  to  300*C,  one  probe  good  for  measurements  up  to  600*C  and  one 
high-temperature  probe  good  for  measurements  to  1200*C.  However,  the 
high  temperature  probe  received  was  not  satisfactory,  and  the  high  tempera- 
ture probe  used  was  assembled  by  Broadview  Research  Corporation  per- 
sonnel. 

The  measurement  apparatus  for  resistance  thermometers  consisted 
essentially  of  a Wheatstone  bridge  called  a "Thermo-Graf"  and  manufactured 
by  the  Hallikainen  Instrument  Company  of  Berkeley,  California.  The  bridge 
required  for  high-accuracy  measurements  must  have  great  resolving  power 
and  stability.  When  a long  lead  intervenes  between  the  probe  and  the  bridge, 
it  is  necessary  to  compensate  for  the  resistance  of  the  lead.  Moreover,  when 
ac  is  used,  it  is  necessary  also  to  compensate  for  the  capacitance  of  the  leads. 

Several  sensitivity  ranges  were  available;  in  the  most  sensitive, 
precision  was  0.001’C.  The  Thermo-Graf  uses  an  arbitrary  scale  system 
of  decades  whose  interrelation  was  established  experimentally.  Each  of  the 
probes  was  calibrated  against  the  Thermo-Graf  at  four  temperatures.  This 
work  was  done  at  the  laboratories  of  Hallikainen  Instruments,  through  the 
courtesy  of  Mr.  Hallikainen.  A large  paraffin  oil  bath  with  adequate  stirring 
was  used.  The  temperature  standards  were  special  expanded- scale  mercury- 
in-glass  thermometers,  calibrated  against  primary  standards  in  the  labora- 
tories of  the  Shell  Development  Co.,  Emeryville,  California,  or  certified  by 
the  Nation  Bureau  of  Standards. 
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By  the  use  of  a three-lead  circuit,  the  effect  of  the  cable  resistance 
was  eliminated  on  the  assumption  that  variations  in  the  resistances  of  the 
leads  in  the  cable  were  identical.  The  capacitance  of  1,  000  feet  of  cable  was 
compensated  for  by  capacitors  in  the  bridge  circuit  of  the  Thermo-Graf. 

For  less  accurate  measurements  the  resistance  of  the  platinum  ther- 
mometer could  be  determined  to  about  0. 1%  by  using  a portable  dc  Wheat- 
stone bridge  manufactured  by  Leeds  and  Northrup  Company.  Accurate  meas- 
urements (to  * 0. 01*C  or  better)  were  to  be  determined  using  the  ac  bridge 
mentioned  above. 

2.  3.  2 Field  Calibrations  and  Operations 

Certain  exigencies  of  field  types  of  operations  raised  the  need  for 
recalibrations  in  the  field  and  modification  of  original  plan  of  attack.  These 
recalibrations  and  changes  are  described  in  the  following  sections. 

No  further  calibration  was  done  with  the  ac  bridge,  but  probe  failures 
forced  the  calibration  of  additional  equipment  for  use  with  the  dc  bridge.  This 
was  done  in  the  field  at  the  Nevada  Test  Site  by  means  of  adapting  a dewar 
flask  for  use  as  a comparator  bath.  A BKH  #61610  mercury  thermometer, 
having  a range  -20*  to  110*C  with  an  error  of  ± 0.  3*C  for  the  entire  range, 
was  used  to  measure  the  temperature  of  the  bath.  The  dewar  was  fitted 
with  a 750-watt  heating  element  and  electric  stirring  device.  A block  of 
styrofoam  insulation  was  added  as  a lid  to  keep  thermal  losses  at  a minimum. 
The  bath  was  filled  with  distilled  water  which  was  heated  in  15*  steps  to  approx 
imately  90*C.  With  this  apparatus,  it  was  possible  to  maintain  the  tempera- 
ture of  the  bath  to  * 0. 1*C  while  the  platinum  resistance  was  measured.  All 
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probes  were  calibrated  with  300  or  1000  feet  of  three-conductor  copper 
cable  attached. 

The  field  calibrated  probes  when  cross  calibrated  with  the  original 
Broadview  Research  probes  were  found  to  be  within  ± 0.  6%  of  the  Broad- 
view Research  results.  The  reliability  of  this  field  calibration  varied  from 
± 0.  4*  at  0*C  to  ± 0.  5*  at  90*C. 

After  the  measurements  were  completed,  a temperature  bath  was 
constructed  to  recalibrate  the  probes.  The  reliability  of  this  calibration 
varied  from  1 0.  26*  at  0*C  to  * 0.  38*  at  90*C  and  agreed,  as  shown  in 
Table  2.  1 with  the  earlier  field  calibrations. 


TABLE  2.  1 COMPARISON  OF  PLATINUM  RESISTANCE 
THERMOMETER  CALIBRATIONS 


Resistance  of 
platinum,  ohms 

Temperature,  *C 
field  calibration 

Temperature,  *C 
laboratory  calibration 

Temperature 
difference,  *C 

109.  37 

22.69 

22.58 

-0. 11 

112.45 

30.  54 

30.  38 

- .16 

116.48 

40.81 

40.67 

- .14 

120. 11 

50.06 

49.87 

- .19 

124. 13 

60.  31 

60. 17 

- .14 

127.87 

69.84 

69.67 

- .17 

131.93 

80. 19 

80. 12 

- .07 

135.55 

89.42 

89.  37 

- .05 

This  bath  consisted  of  a 5.  5 gallon  container  of  distilled  water  which 
was  indirectly  heated  by  eight  200-watt  light  globes  equally  spaced  around 
the  container  and  supported  within  a thermally  insulated  box  3 ft  square. 
Temperature  control  was  accomplished  by  immersing  in  the  liquid  a ther- 
mistor connected  to  one  arm  of  an  ordinary  Wheatstone  bridge.  The  re- 
sistance variation  of  the  thermistor  with  temperature  determined  a bridge 
voltage  output.  This  output  was  then  fed  to  an  amplifier  which  controlled 
a motor-driven  rheostat  in  the  bath  heater  circuit.  The  temperature  was 
controlled  to  better  than  ± 0.07*C  with  this  apparatus. 

To  measure  the  bath  temperature,  a Fisher  #N1135  mercury  ther- 
mometer calibrated  by  the  Bureau  of  Standards  was  used.  This  thermometer 
had  a range  of  -5*  to  101*C  with  a maximum  error  of  ± 0.07*C  over  the 
entire  range. 

For  diagnostic  purposes,  four  holes  were  drilled  into  the  shot  area. 
Table  2.  2 and  Figs.  2.4,  2.5,  and  2.6  describe  these  holes. 

Temperature  measurements  were  made  in  these  holes  with  the  platinum 
resistance  probe  and  the  D.C.  Wheatstone  Bridge  or  the  ac  Resistance  bridge. 
Referring  to  Fig.  2.7,  resistance  values  for  R|,  R£,  R-,  R.„,  and  R__  had 
been  determined  at  a reference  temperature  Tq  where: 

Rj,  R2  and  R ^ were  resistances  of  the  copper  conductors. 

Rc  was  the  resistance  of  a manganin  resistor. 

Rpt  was  the  resistance  of  the  platinum  resistor. 


25 


TABLE  2.  2 DESCRIPTION  OF  LOGGING  HOLES 


Hole 

No. 

Drilling  Station 

Total 
depth,  ft 

Total  ft  of 
casing  in  hole 

Closest  distance 
from  hole  to  G.  Z. 

1 

Mesa 

912* 

897  NXT 

7 ft  above.  19  ft  N 32*W 

B 

In  tunnel,  211.6ft 
from  G.  Z. 

250.5 

250.  7 NX 

56.  8 ft  below,  12  ft  right 

C 

In  tunnel,  211.6  ft 
from  G.  Z. 

284.5 

284.  5 NX 

27.  8 ft  below,  26  ft  right 

D 

— * — 

In  tunnel,  211.6  ft 
from  G.  Z. 

288.5 

211.5  NX+t 
288.  5 BXTT 

3 ft  above,  4 ft  right 

lodged  in  the  bottom  of  the  hole. 

^ 3 inch  i.  d.  x l/4  in.  wall  thickness. 

^ 2-3/4  in.  i.  d.  x 3/l6  in.  wall  thickness. 


lABx  = Rc  + R1t 
o Ao 


+ R, 


lBCT  = RPtT  + R2_  + R3_ 
Ao  Ao  o Ao 


(1) 

(2) 


The  manganin  resistor  potted  in  the  probe  had  a very  low  temperature  coeffi- 
cient of  resistivity.  Therefore!  at  the  temperatures  encountered  (0*  to  100*C) 
it  had  a constant  resistance  R . To  determine  Rp  (see  Fig.  2.  7): 


RPt  = RBC„ 


- (R,  +R,  ) 


lPt  = RBCt  " |*R2_  + R3_  * + A(R2  + R3)T 

T L To  To  °J 


RPt  = RBC„ 
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Since  all  three  copper  conductors  lie  in  the  same  cable  and,  hence,  had  the 
same  temperature  distribution,  the  percentage  change  of  resistance  from  T 
to  T would  be  equal  in  Rp  R^,  and  R y Thus:  c 


a(r2  + r3)J  a(r1  + r2)J 

o _ o 
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Therefore 


RPt  = rbct  " 

All  of  the  parts  of  the  right  hand  side  of  equation  (7)  could  be  measured 
directly. 

The  ac  resistance  bridge  was  used  to  measure  temperatures  up  to  25 *C. 
Higher  temperatures  could  not  be  measured  because  the  bridge  calibration 
was  upset  by  capacitance  variations  with  temperature  in  the  300-foot  3-con- 
ductor cable.  As  seen  in  Fig.  2.  7.  the  probe  connects  to  the  ac  bridge  such 
that  the  cable  lead  resistance  external  to  the  probe  appears  in  opposite  arms 
of  the  bridge.  Thus,  the  platinum  resistance  may  be  determined  by  using 
equation  (8)  derived  from  the  Broadview  Research  calibration. 

•C*  0.  0055287  (total  divisions)  - 19.  489  + Pc  (8) 

where  P = individual  probe  correction.  The  values  for  total  division  could 
be  read  Sirectly  from  the  bridge. 

The  platinum  resistance  probe  was  positioned  in  the  holes  using  1/8- 
in.  steel  pipe  that  was  thermally  insulated  from  the  sensing  element,  and  the 
time  required  for  the  probe  to  reach  thermal  equilibrium  was  determined. 

Since  it  required  up  to  one  to  two  hours  for  the  probes  to  reach  equilibrium, 
two  methods  were  used  in  making  the  temperature  measurements.  During 
the  drilling  process,  10-minute  measurements  were  made  at  the  bottom  of 
the  hole  as  the  hole  was  advanced  and  converted  to  equilibrium  temperatures. 
This  was  done  to  delay  the  drilling  as  little  as  possible  but  yet  obtain  measure- 
ments as  the  hole  advanced.  After  the  holes  were  complete,  the  probe  was 
positioned  in  the  holes  at  5-,  10-,  or  20-ft  intervals  and  left  at  each  position 
until  equilibrium  was  reached. 

To  obtain  the  normal  tuff  temperatures,  a number  of  holes  were  drilled 
to  depths  of  6 to  40  ft  in  the  walls  of  the  tunnel  at  approximately  200-ft  in- 
tervals from  the  drilling  room  to  the  portal  (see  Fig.  2.  6).  The  probes  were 
sealed  in  these  holes  and  measurements  were  made  two  or  three  times  a 
week  with  the  ac  resistance  bridge. 

Measurements  were  also  made  in  the  vertical  mesa  hole  with  thermis- 
tors by  USGS. 
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CHAPTER  3 


RESULTS 


Rainier  shot  of  Operation  Plumbob  was  fired  at  09  hours.  59  minutes, 
59.45  seconds  on  September  19,  1957  (Reference  4).  The  measurements 
recorded  during  and  after  the  detonation  in  the  installations  described  in 
Chapter  2 are  discussed  in  the  following  sections. 

3.  1 PRESSURES  BEHIND  TUNNEL  BLAST  DOORS 

The  shot  tunnel  was  reentered  within  a few  hours  after  detonation.  The 
pressure  gage  stations  at  the  mouth  of  the  tunnel  indicated  negligible  pressure 
rises  behind  each  of  the  blast  doors  in  the  tunnels.  Upon  penetrating  the  tun- 
nel about  four  hours  after  detonation,  the  blast  doors  were  found  to  be  satis- 
factorily sealed  and  they  pressure  rises  behind  them,  as  read  at  the  local  gage 
stations,  was  below  the  sensitivity  of  the  gages  to  record.  These  observa- 
tions together  with  the  lack  of  radioactivity  in  the  tunnel  indicated  that  the 
tunnel  had  sealed  satisfactorily,  according  to  predictions. 

3.  2 TRANSIENT  TEMPERATURE  RISE  IN  TUFF  DURING  PASSAGE  OF 
SHOCK  WAVE 

Five  probes  were  installed  at  distances  of  approximately  25,  45.  5,  95.  5, 
400,  and  600  ft  from  the  detonation  point  and  grouted  in  place  with  a special 
grout  supplied  by  UCRL.  The  full-scale  deflections  chosen  for  each  of  these 
locations  were  100,  25,  0.5,  0.2,  and  0.  2*C,  respectively. 

Tbe  first  signal  whiCR  appeared  on  the  records  was  essentially  coinci- 
dent with  the  detonation  of  the  device  and  apparently  consisted  of  an  electro- 
magnetic disturbance  which  was  initiated  by  the  detonation.  On  the  more 
sensitive  ranges  (the  two  most  distant  sensors)  the  signal  was  of  such  a mag- 
nitude that  the  recording  apparatus  shows  a permanent  "set"  so  that  the  signal 
saturated  the  amplifiers  for  a long  period  of  time  and  obscured  any  tempera- 
ture signals  which  might  have  been  present.  Although  the  existence  of  this 
spurious  signal  had  been  anticipated  and  the  gages  were  as  well-shielded  as 
was  felt  compatible  with  the  desired  temperature  measurement,  the  signal 
was  large  enough  to  completely  cover  and  mask  the  desired  information  on 
these  channels. 

The  three  close-in  channels  (25,  45.5,  and  95.5  ft,  respectively)  were 
recovered  partially  after  the  spurious  signal,  but  the  signals  which  followed 
were  not  interpretable  in  terms  of  temperature.  One  of  the  reasons  for  this 
is  the  low-frequency  response  of  the  system  which  did  not  permit  resolution 
of  information  containing  frequency  components  higher  than  about  500  cycles 
per  second.  It  is  possible  also  that  strain  effects  may  have  masked  the  de- 
sired temperature  signal  even  though  precautions  to  exclude  these  effects 
were  taken  in  the  design  of  the  temperature  sensors.  In  addition,  low-fre- 
quency electromagnetic  signals  may  have  persisted  long  enough  to  cause 
confusion  in  the  recorded  signals.  The  signal  recorded  on  one  of  the  three 
close-in  channels  is  reproduced  in  Fig.  3. 1. 

3.  3 POST -SHOT  TEMPERATURE  MEASUREMENTS  IN  THE  MEDIUM  SUR- 
ROUNDING THE  NUCLEAR  DETONATION 

The  temperature  distributions  in  the  four  diagnostic  holes  are  shown 
in  Figs.  3.  2,  3.  3 and  3. 4.  A plot  of  the  isothermal  lines  on  a surface  de- 
fined by  the  holes  B,  C,  and  D is  shown  in  Fig.  3. 5.  The  highest  tempera- 
ture observed  was  89.  7 * 0.4*C.  The  normal  tuff  temperature  at  this  depth 


- 35 


/ 


4 


. frvv 

J " 


L 


is  16.  6*0.  1*C.  Heating  from  the  shot  occurred  as  far  as  180  ft  from  ground 
zero  (G.Z.)  in  the  B,  O,  and  D holes  and  extends  approximately  87  ft  above 
G.Z.  in  the  vertical  .'  ole.  This  asymmetry  is  clearly  shown  in  Fig.  3.5.  This 
figure  indicates  that  t-ie  center  of  the  heat  source  is  located  approximately  50 
ft  below  G.Z.  and  that  the  maximum  temperature  is  less  than  or  equal  to  the 
boiling  point  of  water  (approximately  94*C  at  shot  altitude). 

The  exact  boundary  of  the  heated  zone  is  difficult  to  determine  because 
of  the  thermal  perturbations  introduced  by  the  drilling  operation.  These  per- 
turbations include  the  heating  of  the  rock  by  the  drill  bit,  the  cooling  by  the 
drill  water  and  the  heat  transfer  in  the  casing  in  the  hole.  These  effects  are 
shown  in  Fig.  3. 5 in  the  region  below  20*C. 

Temperature  measurements  were  made  at  three  points  in  a vertical  hole 
drilled  from  a point  on  the  mesa  899  ft  above  and  1 3.  5 ft  N 48*  30'  W of  G.  Z. 

At  539  ft  the  temperature  measurement  was  12.  2 dk  0.4*C,  at  670  ft  it  was 
11.  7 * 0.4*C  and  at  685  ft  it  was  10.  2 * 0.4*C.  The  temperature  at  539  ft 
was  measured  at  the  bottom  of  an  inverted  cone-shaped  cavity  26  ft  high  and 
50  ft  in  diameter.  This  cavity  resulted  from  the  collapse  of  the  initial  cavity 
formed  at  G.  Z.  The  12.  2 ± 0. 4*C  temperature  measured  at  this  point  com- 
pares with  a pre-shot  temperature  of  12.6  * 0.  1*C  recorded  by  USGS  at  the 
same  depth  in  a vertical  hole  located  87  ft  N 65*20'  W from  this  post-shot  hole. 
The  result!  observed  at  670  ft  and  685  ft  are  lower  than  the  pre-shot  values 
obtained  by  USGS  (Reference  5)  at  the  same  depth  (Figs.  3.6  and  3.  7).  This 
is  believed  to  be  a result  of  the  large  amount  of  drilling  water  lost  to  the  me- 
dium during  the  drilling  process.  Two  temperature  logs  of  this  hole  by  USGS 
(Reference  5)  20  days  apart,  two  months  later,  indicated  that  the  temperature 
at  these  points  was  still  below  but  approaching  pre- shot  values  (see  Fig.  3.7). 

During  the  post-shot  drilling  operation,  approximately  68,000  gallons 
of  water  (10*-15*C)  were  lost  in  the  vertical  mesa  hole  between  513  fit  and  899 
ft  and  approximately  7,000  galloiu  (10*-15*C)  in  holes  B,  C,  and  D.  This 
is  small  compared  with  the  3(10) 7 gallons  (20%  by  weight)  of  water  contained 
in  the  pre-shot  tuff  in  the  heated  zone.  Therefore,  the  gross  effect  in  cooling 
the  heat  source  is  negligible.  However,  the  effect  can  be  seen  in  the  boundary 
areas  as  shown  in  a comparison  of  the  pre-shot  and  post-shot  USGS  measure- 
ments (Figs.  3.6  and  3.  7).  The  water  lost  to  the  medium  in  drilling  the  B, 

C,  and  D holes  was  lost  in  the  region  within  60  ft  of  G.Z.  This  corresponds 
to  the  initial  cavity  formed  by  the  shot. 

The  ambient  rock  temperature  along  the  tunnel  was  measured  in  6-  and 
40-ft  deep  horizontal  holes  as  shown  in  Fig.  3. 8.  A temperature  of  16.  6 * 

0. 1*C  was  measured  at  the  G.Z.  depth.  T'e  temperature  along  the  tunnel 
decreases  as  the  portal  is  approached  as  shown  in  Fig.  3.8.  The  6-ft  hori- 
zontal holes  were  not  satisfactory  for  determining  background  rock  tempera- 
tures because  of  the  cooling  resulting  from  the  air  circulating  in  the  tunnel. 

In  the  40-ft  horizontal  holes,  a cooling  of  0.05  to  0. 1*C  was  observed  in  a 
period  of  56  days. 

After  the  drilling  was  completed,  two  temperature  logs  were  made  in 
the  B hole  58  days  apart  (144  and  202  days  after  the  shot).  As  seen  in  Fig. 

3.  2,  these  curves  have  essentially  the  same  shape.  Integrating  the  area 
under  the  curves,  however,  showed  that  a 6.4%  reduction  had  occurred  in 
58  days.  The  peak  temperature  in  this  period  had  dropped  from  89.  7*C  to 
84.  2*C.  A similar  comparison  of  two  temperature  logs  made  45  days  apart 
(118  and  163  days  after  the  shot)  in  the  D hole  showed  a reduction  of  3.  2% 
in  the  area  of  overlap.  The  average  rate  of  change  in  the  B hole  is  -0. 1 1 * 
0.02%  per  day.  The  average  rate  of  change  in  the  D hole  is  -0.07  * 0.02% 
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per  day*  Figure  3.2  also  shows  a partial  log  made  in  the  B hole  before  casing 
the  hole  and  during  the  drilling  operation.  It  is  difficult  to  interpret  the  com- 
parison of  this  logging  with  that  of  the  cased  hole  because  not  enough  informa- 
tion is  known  on  the  effects  of  casing,  drilling  water,  etc. 

A tunnel  is  now  being  driven  into  a point  approximately  20  ft  vertically 
below  G.Z.  From  this  point,  more  diagnostic  holes  will  be  drilled  to  ob- 
tain better  definition  of  the  geometry  and  magnitude  of  the  shot  effects. 


i 


HOLE  C— RAINIER  TUNNEL 


Figure  SJ 


>001 


DEPTH  , feet 


Figure  3.6 


400 


500 


600 


SURFACE 


TEMPERATURE,  °C 
10  12  14  16 


CABLE  321 
4 8/9/57 

• 8/8/57 

* 8/13/57 


DEPTH  t feet 


°°oo ooo  o 8 g 


Figure  3.8 


CHAPTER  4 


DISCUSSIONS  AND  RECOMMENDATIONS 


The  results  of  the  measurements  described  in  the  previous  chapters, 
when  viewed  together  with  results  obtainable  from  the  other  experiments 
conducted  on  this  shot,  give  a coherent  picture  of  the  sequence  of  events  in 
the  nuclear  detonation  fired  deep  underground  in  Nevada  tuff.  The  general 
conclusions  which  can  be  drawn  from  the  results  of  the  measurements  des- 
cribed in  the  previous  chapter  are  discussed  below. 


4. 1 EFFECTIVENESS  OF  TUNNEL  DESIGN 

The  lack  of  static  pressure  in  the  tunnels  behind  the  blast  doors  and 
also  the  lack  of  damage  to  the  blast  doors  indicated  that  the  original  concept 
of  the  tunnel  design  was  effective  in  preventing  the  energy  of  the  detonation 
or  explosive  products  from  venting  into  the  tunnel. 

It  is  recommended  that  this  blast-door  instrumentation  be  included  in 
any  similar  type  futnre  installations. 

4.  2 TEMPERATURE  RISE  IN  TUFF  BEHIND  SHOCK  FRONT 

No  direct  results  were  obtained  on  this  quantity,  as  discussed  in 
Section  3.  2 of  Chapter  3.  However,  the  long-term  post-shot  temperature 
measurements  indicate  that  most  of  the  heat  energy  is  confined  to  a region 
less  180  ft  in  radius  from  the  original  point  of  detonation.  This  meas- 
urement indicates  that  the  general  concept  of  rapid  attenuation  of  pressure 
with  radius  and  the  accompanying  reduction  in  temperature  rise  in  the  shock 
wave  is  valid.  As  was  seen  in  Section  1.  3.  2 no  appreciable  temperature 
rise  was  anticipated  for  shock  pressures  of  10^  bars  or  less,  correspond- 
ing to  radii  of  less  than  100  ft  from  the  explosion  point. 

The  measurement  should  probably  only  be  repeated  when  some  ad- 
vance in  instrumentation  technique  renders  the  success  of  the  measurement 
more  probable. 

4.  3 POST-SHOT  TEMPERATURE  MEASUREMENTS  IN  THE  MEDIUM  SUR- 
ROUNDING THE  POINT  OF  DETONATION 

The  measurements  indicate  that  no  temperatures  greater  than  the  boil- 
ing point  of  water  at  these  elevations  (94*C)  exist  anywhere  in  the  region  sur- 
rounding the  detonation  point.  The  hot  material  is  collected  in  a dish- shaped 
region  about  55  to  60  ft  below  the  original  detonation  point.  Interior  to  and 
above  this  material  is  relatively  unconsolidated  soil  and  rock,  whereas  below 
the  hot  region  is  a layer  of  crushed  and  deformed  tuff.  The  highest  tempera- 
ture was  reached  in  the  center  of  vitreous  material,  which  appears  to  have 
been  molten  tuff  (see  Fig.  3. 5). 

A rough  integration  over  the  temperature  distributions  given  in  Fig.  3.5 
may  be  made  as  follows:  Assume  that  the  temperature  distributions  are 
roughly  symmetrical  about  an  axis  passing  through  the  original  ground  sero 
and  the  highest  temperature  region  encountered.  The  volumes  inclosed  by 
the  constant  temperature  contours  may  then  be  approximated  by  hemispheres 
having  radii  as  estimated  in  Table  4. 1.  Using  an  average  specific  heat  of 
0.  35  calories  per  gram,  a density  of  1.9  grams  per  cc,  we  calculated  as  in 
Table  4. 1 that  the  total  amount  of  energy  contained  in  the  tuff  in  the  form  of 
material  heated  above  20*C  is  8. 6 X 10“  calories,  or  over  half  of  the  total 
bomb  energy  of  1.  7 X lO1^  calories. 
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TABLE  4. 1 CALCULATION  OF  TOTAL  HEAT  CONTENT  OF  TUFF 
Temp.,  *C  Radius,  cm  Vol. -1^  Temp,  rise,  Heat  content  «p  VCp  At. 


20*C 

3000 

5.67  X 1010 

3. 4*C 

12.8  X 10 

30  *C 

2670 

3.98 

10.0* 

26.1 

40*C 

2460 

3. 12 

10. 0» 

20.8 

50'C 

2160 

2.12 

10.0* 

14.1 

60*C 

1770 

1.16 

10.0* 

7.7 

70*C 

1320 

0.483 

10.0* 

3.2 

80*C 

1020 

0.  222 

10.0* 

1.5 

86. 2 X 10J 


An  explanation  for  the  above  measurements  which  is  consistent  with 
all  the  data  is  as  follows: 

Shortly  after  detonation  a molten  or  fused  shell  was  formed  on  the  in- 
side of  the  shot  chamber.  Almost  simultaneously  the  walls  of  the  chamber 
underwent  hydrodynamic  motion  until  a cavity  with  dimensions  of  50  to  60  ft 
was  formed.  The  amount  of  material  in  the  fused  layer  may  be  estimated 
from  radiochemical  measurements,  since  essentially  all  of  the  fission  frag- 
ments were  trapped  in  the  vitrified  layer  at  the  edge  of  the  cavity.  The  vol- 
ume of  fused  material  (fusing  temperature  1 200-1 500 *C)  is  estimated  to  be 
3. 5 x 10®  cm®,  corresponding  to  a shell  thickness  of  about  4 in.  at  the  cavity 
radius  of  55ft.  Shortly  after  the  detonation  the  roof  of  this  cavity  collapsed 
in  such  a manner  that  most  of  the  roof  material  fell  to  the  bottom  of  the  cav- 
ity and  the  cavity  filled  with  unconsolidated  permeable  material.  This  col- 
lapse was  propagated  upward  until  a stratum  capable  of  supporting  the  stress 
was  encountered. 

The  hot  material  was  assume  A to  cool  rapidly  to  a temperature  below 
100*C  due  to  the  admixture  of  unconsolidated  tuff,  when  the  cavity  roof  struc- 
ture caved  in,  and  by  heating  of  the  arge  quantity  of  water  contained  in  this 
mass  of  rock  as  pore  water  or  water  of  crystallisation.  The  boiling  point 
for  water  at  this  elevation  is  about  94*C.  Taking  into  account  the  fact  that 
the  once  molten  material  is  thee  mixed  essentially  with  a few  feet  of  radial 
thickness  of  unconsolidated  tuff,  and  adding  a small  amount  of  radial  tuff 


thickness  to  account  for  material  shock  heated  to  shove  94* C,  the  total  mass 
of  hot  tuff  at  a 55-ft  radius  is  seen  to  be  about  10*°  grams,  containing  about 
2 x 10?  g of  water.  To  evaporate  this  amount  of  water  would  require  2 X 10“  X 
[100-16  + 540]  « 1. 24  X 10*2  cal.  The  total  bomb  energy  is  approximately 
1. 7 X 1012  cal  so  that  is  is  quite  reasonable  to  assume  that  any  high  tempera- 
tures in  the  tuff  were  soon  reduced  to  below  the  region  of  90  to  100*C.  In  a 
similar  fashion,  by  assuming  400  cal/g  as  the  latent  heat  of  fusion  for  tuff 
and  the  estimate  of  7 X 10®  g of  fused  material  for  radiochemical  measure- 
ments, the  amount  of  heat  which  entered  the  fused  material  is  estimated  as 
about  30%  of  the  total  bomb  energy. 
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The  feet  that  the  broken  unconsolidated  material- is  permeable  to  water 
tends  to  invalidate  the  procedural  cc  icept  of  long-term  temperature  measure- 
ments of  heat  conduction  in  the  tuff,  since  heat  will  be  transferred  by  water, 
vapor  or  even  gas  conduction  and  convection  to  regions  of  the  tuff.  For  this 
reason  very  accurate  (0. 001°C)  long-term  measurements  of  the  heat  diffusion 
by  conduction  away  from  the  central  source  will  not  yield  easily  interpretable 
results  and  should  probably  be  discontinued.  On  the  other  hand,  it  is  vital 
to  document  the  sequence  of  events  in  the  shocked  medium  by  moderately  ac- 
curate temperature  logging  of  drill  holes  made  before  and  after  the  nuclear 
detonation. 


APPENDIX  A 

This  appendix  gives  the  solution  of  the  integral 

f*  r.  . .-<»r>2/4.«  it, 
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Similarly,  the  second  integral  is  solved  by 
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The  integral  becomes  2st  £ .-Ydr . ,r  .-(^')2/4.«  dr, 
and  its  solution  is 
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Combining  the  two  we  obtain 
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